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v«D . Abstract 

(N ■ 

^ ■ 

Q ' SN 2002ic is a unique supernova which shows the typical spectral features 

^ ■ of Type la supernovae (SNe la) near maximum light, but also apparent hydrogen 

"^ ' features that have been absent in SNe la. We have calculated hydrodynamical 

Qh' models for the interaction between the SN la ejecta and the H-rich circumstellar 

O ' medium (GSM) to reproduce the observed features of SN 2002ic. Based on our 

-^ ' modeling, we suggest that GSM is aspherical (or highly clumpy) and contains ~ 

^ ■ 4-5 Mq. Possible progenitor systems of SN 2002ic are discussed. 

>'■ 

K> , 1. Type la Supernovae and Circumstellar Medium 

^ : 

■ - - ' Type la supernovae (SNe la) are characterized by the lack of hydrogen and 

the prominent Si line in their spectra near maximum light and widely believed 
to be thermonuclear explosions of mass-accreting white dwarfs in binary systems. 
SNe la have been used as a "standard candle" to determine cosmological param- 
eters thanks to their relatively uniform light curves and spectral evolution. SNe 
la are also the major sources of Fe in the galactic and cosmic chemical evolution. 
Despite such importance, the immediate progenitor binary systems have not been 
clearly identified yet (e.g., [TIHIT!^). 

For a model of SN la progenitors, Hachisu et al. E] proposed a single 
degenerate model in which the white dwarf blows a massive and fast wind (up 
to 10~^Mq yr~^ and 2000 km s"^) and avoids a formation of common envelope 
when the mass transfer rate from the normal companion exceeds a critical rate of 
~ 1 X 10~® Mq yr~^ pi]. Such an evolutionary phase is dubbed "accretion wind 
evolution" instead of "common envelope evolution." Such a binary can keep its 
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separation almost unchanged. The white dwarf can steadily accrete a part of the 
transferred matter and eventually reach the Chandrasekhar mass. 

In the strong wind model, the white winds form a circumstellar envelope 
around the binary systems prior to the explosion. When the ejecta collide with 
the circumstellar envelope, X-rays, radio, and Ha lines are expected to be emitted 
by shock heating. Attempts have been made to detect such emissions, but so far 
no signature of circumstellar matter has been detected. 

The upper limit set by X-ray observations of SN 1992A is M/viq = (2 — 3) 
X 10^^ Mq yr~^ pn]- Radio observations of SN 1986G have provided the most 
stringent upper limit to the circumstellar density as M/viq = 1 x 10"^Mq yr~^, 
where fio means fio = w/10 km s~^. This is still 10 — 100 times higher than 
the density predicted for the white dwarf winds, because the white dwarf wind 
velocity is as fast as ~ 1000 km s~^. For Ha emissions, the upper limit of M/viq = 
6 X 10-6 Mq yr-i has been obtained for SN 1994D. 

2. SN 2002ic 

SN 2002ic was discovered on 2002 November 13 UT at magnitude 18.5 by 
the Nearby Supernova Factory search [THj. Hamuy et al. jHj reported strong Fe 
III features and a Si II A6355 line in the early-time spectra of SN 2002ic and 
classified it as a SN la. 

However, strong Ha emission was also observed. The emission was broad 
(FWHM > 1000 km s^^) suggesting that it was intrinsic not to an H II region of 
the host galaxy but to the supernova. The detection of Ha is unprecedented in a 
SNIa(e.g., [imn]). 

Hamuy et al. |H] suggested that it arose from the interaction between the 
SN ejecta and a dense, H-rich circumstellar medium (GSM), as in Type Iln SNe 
(SNe Iln). If this interpretation is correct, SN 2002ic may be the first SN la 
to show direct evidence of the circumstellar (CS) gas ejected by the progenitor 
system, presenting us with a unique opportunity to explore the GSM around a 
SN la and the nature of the progenitor system. 

2.1. Spectroscopic Features of SN 2002ic 

The late-time spectrum of SN 2002ic is strikingly similar to those of Type 
Iln SNe 1997cy ^7] and 1999E ^^ as shown in Figure 1. Spectroscopically, 
SN 2002ic evolve with time significantly, in particular in the Ha line and its 
complex profile. Hamuy et al. jH] detected in the early-time spectra an unresolved 
narrow Ha emission on top of a ~ 2000 km/s base, which were superimposed on 
dominant SN la line features. One year after the explosion, however, the Ha line 
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Fig. 1. Spectral comparison between SNe 2002ic (red thick lines; Subaru, ~ 222 d), 
1997cy (blue thin lines) and 1999E (green dashed line). 



became much more prominent and consisted of a narrow core and a ~ 5000 km/s 
component. Other strong features identified in Figure 1 include Ca and O lines 
as broad as > 10, 000 km/s and broad permitted Fe II multiplets [3]. 

SNe 1997cy and 1999E were initially classified as Type Iln because they 
showed Ha emission. SN 2002ic would also have been so classified, had it not been 
discovered at an early epoch. SN 1997cy [z = 0.063) is among the most luminous 
SNe discovered so far (My < —20.1 about maximum light), and SN 1999E is also 
bright {Mv < -19.5). Both SNe 1997cy and 1999E have been suspected to be 
spatially and temporally related to a GRB [H^^. However, both the classification 
and the associations with a GRB must now be seen as highly questionable in view 
of the fact that their replica, SN 2002ic, appears to have been a genuine SN la at 
an earlier phase. 



2.2. Observed Light Curve of SN 2002ic 

The VVOIR bolometric light curve of SN 2002ic has been constructed by 
Deng et al. [3] from the available BVRI photometry and the spectrophotom- 
etry |H1 E] as shown in Figure 2. To construct the light curve of SN 2002ic, 
we first integrated the Subaru spectrum. This yielded L = (5.9 ± 0.6) x 10^^ 
ergs s~^, corresponding to Mboi ~ —18.2, assuming a distance of 307 Mpc {Hq = 
65 km s~^ Mpc~^). The bolometric corrections thus estimated was used to con- 
vert the early-time photometry in Hamuy et al. jS] and the late-time MAGNUM 
telescope photometry into bolometric luminosities. 



The light curve of SN 2002ic is brighter at maximum and dechnes much 
more slowly than typical SNe la jH]. The late time light curve of most SNe is 
powered by the radioactive decay of ^^Co to ^^Fe. The decline of SN 2002ic is 
much slower than the Co decay rate, which indicates the presence of another 
source of energy. 

In fact the overall light curve of SN 2002ic resembles SNe 1999E ^3] and 
1997cy HZI (see Figure 2). We use UBVRI bolometric light curves of SNe 1997cy 
and 1999E for comparison [T7lll4j . Assuming E{B — V) = 0.06 for the Galactic 
extinction (NED), SN 2002ic is only a factor of 1.3 dimmer than SN 1997cy, but 
more than 100 times brighter at late phases than typical SNe la. The light curve 
of SN 1997cy has been modeled in the context of circumstellar interaction [T7] . 
which is very likely the same energy source for SN 2002ic. 



43.5 



43.0 






42.0 



41.5 



41.0 



■ 1^ 


' ' ' 


' 


' ' ' 


' 





. ( ' 

-; ' 

"i ^ 
"' \ 


2©o 


i 


X 


■ 

O 
X 


SN 2002ic 
SN1997cy ■ 
SN1999E 
- SN1994D 

X 


\ 

\ 










X O- 


\ 










X Q 


" .... 1 .... 1 


. . . 


1 1 1 


. . . 


1 , 


, |- 



100 200 300 400 

Rest-frame days 



500 



600 



Fig. 2. Comparison of the UBVRI bolometric light curves of SN 2002ic (red filled 
squares) with those of SNe 1997cy (blue open circles) and 1999E (green crosses), 
and the normal SN la 1994D (black dashed line). 



3. Circumstellar Interaction Models 

We calculated the interaction between the expanding ejecta and CSM (de- 
tails will be seen in Suzuki et al., in preparation). For the supernova ejecta, we 
used the the carbon deflagration model W7 [12]; its kinetic energy is E = 1.3 x 
]^g5i g]^g Pqj, CSM -vve assumed the power-law density distribution: 



p = poir/Ro) "g cm" 



(1) 
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Fig. 3. Model light curve (black thick line) compared with the observation of SN 
2002ic (red filled circles 0). 



where the parameters are the radius (i?o) and density (po) of the point where the 
ejecta and CSM start interacting, and the index (n) of the density distribution of 
CSM. These parameters are constrained from comparison with the observed hght 
curve. The spherical Lagrange hydrodynamical code and input physics are the 
same as in Suzuki & Nomoto [T^ . 

When the expanding ejecta interacts with CSM, the interaction creates 
the forward shock which is propagating through the CSM and the reverse shock 
which is propagating through the ejecta. 

Shocked matter is heated to T ~ 10^ K for the reverse shock and T ~ 
10^ K for the forward shock. Both shocked regions emit thermal X-rays. For 
the reverse shock, because of relatively high densities in the ejecta, cooling time 
scale is shorter than shock propagation so that the shocked ejecta soon forms 
a dense cool shell ^Hj- This dense cool shell absorbs the X-ray and re-emits in 
UV-optical. This re-emitted photons are observed. We assume that a half of the 
X-rays emitted in the reverse-shocked ejecta is lost into the supernova center, and 
that the other half is transferred outwardly through the cooling shell. We also 
assume that a half of the X-rays emitted in the CSM is transferred inwardly to 
be absorbed by the cooling shell. We take into account the change in time of the 
column density of the cooling shell to evaluate the X-rays absorbed by the shell 
and the optical luminosity. 

Figure 3 shows the successful model for the light curve of SN 2002ic. Here 
Ro = 2 X 10^^ cm, Po = 4 X lO^^^ g cm'^, and n = 1.8 for inner CSM of 4 Mq. 
In the early phase, the model with n = 2.0 (steady mass loss) declines too fast to 
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Fig. 4. Velocity profiles in the interacting ejecta and CSM. 



be compatible with the observation. This imphes that CSM around the SN was 
created by unsteady mass loss of the progenitor system. 

After day ~ 350, the hght curve starts declining. To reproduce the declin- 
ing part of the light curve, we add the outer CSM of 0.7 Mq where the density 
declines sharply as n = 6. This implies that the total mass of CSM is ~ 4.7 Mq. 

We note in Figure 4 that the velocity of the ejecta decelerated by the CSM- 
interaction is < 4000 km s~^ and too low for the value observed in the broad 
spectral features (~ 10,000 km s~^). On the other hand, in order to produce high 
enough luminosity to explain the light curve, such a strong interaction between 
the ejecta and CSM should occur. 

To reproduce both the light curve and the observed velocity of SN 2002ic, 
CSM needs to be aspherical. Suppose the CSM is aspherical consisting of a dense 
region and a thin region. The expanding ejecta interacting strongly with the 
dense region can produce high enough luminosity to explain the light curve. On 
the other hand, the ejecta interacting with the thin region can expand still fast 
enough to be consistent with the observed velocities (see also Deng et al. jB]). 
A pre-existing clumpy disk was also suggested by Wang et al. jTH], based on 
spectropolarimetry. 



4. Discussion 



There are two possible progenitor scenarios for SN 2002ic. One is the the 
explosion of the C+0 core of the massive AGB star (SN I-I-1/2), where the wind 
from the AGB star formed the CSM. The other is the explosion of the white dwarf 
in a close binary blowing wind to create the dense CSM (e.g., [Hill])- 



4-1- Type 1+1/2 Supernovae in AGB Stars 

Single star scenario is the explosion of the massive AGB star whose C+0 
core becomes very close to the Chandrasekhar mass. Before explosion, mass loss 
(super-wind) from the star creates a dense GSM. The G+0 core explodes, which 
is called Type 1+1/2 supernova, and interacts with GSM. 

To make this scenario possible, the metallicity of the system should be 
low because low mass loss rate is necessary for the G+0 to grow to reach the 
Ghandrasekhar mass before the envelope is completely lost. Under the solar 
metallicity, SN 1+1/2 have never been observed. Therefore, we can explain the 
rarity of SN 2002ic-like event assuming that only narrow mass range of AGB stars 
can explode as SNe in low metal environment. 

Aspherical GSM is not unexpected for stars approaching the end of the 
AGB. 

4-2. White Dwarf Winds 

Binary star scenario is the explosion of the accreting G+0 white dwarf 
(same as typical SNe la). However, the companion star is massive and the white 
dwarf blows large amount of accreting gas as accretion wind to create the dense 
GSM. In this scenario, the rarity is can be attributed to the fewness of the com- 
panion star massive enough to produce the quite massive GSM. 

As a progenitor of SN 2002ic, we need a GSM of ~ 4.7 Mq. Such a massive 
GSM is possible only when the donor is as massive as 6 — 7 Mq. For the model 
consisting of a white dwarf and a main-sequence companion jH] , the mass transfer 
rate from such a massive main-sequence companion reaches ~ 1 x 10~^Mq yr~^. 
Then the white dwarf blows a wind of ~ 1 x 10~^Mq yr~^ and the mass stripping 
rate becomes several times larger than the white dwarf wind mass loss rate [7j. 

For the symbiotic model consisting of a white dwarf and a red giant or AGB 
star, the wind mass loss rate can also reach ~ 1 x IQ^^Mq yr"^. In symbiotic 
stars, the mass capture efficiency by the white dwarf is observationally estimated 
to be as small as one or a few percent. Therefore, only when a large part of the 
red giant wind or AGB super-wind is captured by the white dwarf, the white 
dwarf can blow a very massive wind of ~ 1 x 1O~^M0 yr~^ or more. Then, the 
mass stripping rate from the red giant or AGB star also reaches several times 
1O~^M0 yr~^ or more. 

Examples of the accretion wind evolution are identified as transient super- 
soft X-ray sources, i.e., the LMG supersoft X-ray source RX J0513.9— 6951 and 
its Galactic counterpart V Sge [Zj. Especially in V Sge, a very massive wind of 
~ 1 X 1O~^M0 yr~^ has been observationally suggested by the detection of radio. 



Furthermore, the white dwarf wind colhdes with the companion and strips heavily 
off its surface matter. This stripping rate reaches a few or several times the wind 
mass loss rate of the white dwarf, i.e., ~ 1 x 1O~^M0 yr~^ or more [7j. The matter 
stripped off has a much lower velocity than the white dwarf wind itself and forms 
an excretion disk around the binary. Thus the model predict the coexistence of a 
fast white dwarf wind blowing mainly in the pole direction and a massive disk or 
a torus around the binary. Deng et al. jH] propose a new classification. Type Ila 
supernovae, for these events. 
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